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Abstract 

We report on a scheme which might make it practically possible to engineer the 
effective competing nonlinearities that on average govern the light propagation in 
quasi-phase-matching (QPM) gratings. Modulation of the QPM period with a sec- 
ond longer period, introduces an extra degree of freedom, which can be used to 
engineer the effective quadratic and induced cubic nonlinearity. However, in con- 
trast to former work here we use a simple phase-reversal grating for the modulation, 
which is practically realizable and has already been fabricated. Furthermore, we de- 
velop the theory for arbitrary relative lengths of the two periods and we consider the 
effect on solitons and the bandwidth for their generation. We derive an expression 
for the bandwidth of multicolor soliton generation in two-period QPM samples and 
we predict and confirm numerically that the bandwidth is broader in the two-period 
QPM sample than in homogeneous structures. 



Quasi-phase- matching (QPM) is a major alternative over conventional phase 
matching in many laser applications based on frequency- conversion processes 
in quadratic nonlinear media (for reviews, see (0 ||)). Besides other practical 
advantages, QPM allows tailoring the nonlinearity of the material to form com- 
plex structures. This opens a range of new possibilities, which have become ex- 
perimentally feasible with the recent progress in poling techniques. For exam- 
ple, engineerable pulse compression in frequency-doubling schemes in synthetic 
QPM gratings has been demonstrated in aperiodically poled lithium niobate 
and potassium niobate (||; f|; Q ||), and transverse QPM gratings have been 
made both for shaping second-harmonic beams and to extend the spectral cov- 
erage of optical parametric oscillators (|7|; H). Bandwidth enhanced parametric 
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interactions can be obtained in modulated-period structures @, multiple non 
linear interactions can be achieved in quasi-periodic schemes (Jl^; |ll|; 12; 13) 



and simultaneous generation of multiple color laser light has been demon- 
strated in QPM crystals doped with active lasing ions flTil) . QPM engineering 
also finds novel important applications beyond pure frequency-conversion de- 
vices, e.g., to generate enhanced cascading phase-shifts (|15D , all-optical diode 
operation (|TED, and multicolor soliton formation (|T7]; [IJ| [TP]; |2TJ; |2JJ; 



Here we show that two-period QPM structures offer new opportunities for 
soliton control. 

Multicolor solitons mediated by second-harmonic generation (SHG) form by 
mutual trapping of the beams at the fundamental frequency (FF) and at the 
second harmonic frequency (SH). Here we focus on bright solitons whose ba- 
sic properties are well known. Whole dynamically stable families exist above 
a certain power threshold for all phase mismatches. At lowest order the effect 
of QPM grating is to average the quadratic nonlinearity (f^). However, taking 
higher order perturbations into account reveals that the corresponding aver- 
aged field equations include effective cubic, Kerr-like terms(|l9|; |23| ; |2~4f) . Such 



terms modify the average properties of CW waves (B5|; |26|) and the soliton fam- 



ilies of the averaged equations ([19| ; |24[ ; |27|) , which can be analyzed as sustained 



by competing quadratic and effective cubic nonlinearities ( p8| ; p9[ |30| ; |3~T| ; p2|) . 

Nevertheless, in standard QPM the strength of the induced averaged cubic 
nonlinearity is proportional to the ratio between the QPM grating period and 
the soliton characteristic length (i.e., the diffraction length in the case of spa- 
tial solitons). At optical wave lengths, the former is typically of the order 
of ten /im, whereas the latter is a few mm. Therefore, the strength of the 
induced cubic nonlinearities is extremely small. The question thus naturally 
arises whether and how QPM engineering can be employed to bring the effec- 
tive cubic nonlinearities to compete with the quadratic in the average fields. 
One solution is to add a strong dc-part to the nonlinear QPM grating, i.e. 
as done in fl3"3"t p4|) . This adds a term to the induced Kerr terms, which is 
proportional to the QPM grating period and the dc-value squared, and thus 
can be large (|24|). Another potentially more versatile technique is to modulate 
the QPM period with a second longer period, as it was shown theoretically 
in (|2"3]). This introduces an extra degree of freedom, which can be used to 
engineer the effective quadratic and induced averaged cubic nonlinearity. 

Here we consider the latter technique. However, in contrast to the rather 
complicated long-period modulation used in ( [23] ) we use here a simple phase- 
reversal grating for the modulation, similar to the two-period QPM sample 
recently employed by Chou and co-workers (|T^) for multiple-channel wave- 
length conversion in the third telecommunication window. Furthermore, we 
develop the theory for arbitrary relative lengths of the two periods and we 
consider the effect of the modulation on multicolor solitons. We expect soliton 
formation to be possible for a variety of phase-reversal periods, in view of the 
fact that multicolor solitons have been shown to be robust against strong per- 



turbations, including quasi-periodic QPM gratings (pi|) and even periodic gain 
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Fig. 1. Two-period QPM grating. 

and loss(|35l). Using the induced averaged cubic nonlinearities we derive an ex- 
pression for the bandwidth of multicolor soliton generation in two-period QPM 
samples and we show that the bandwidth is broader in the two-period QPM 
sample than in homogeneous quadratic nonlinear materials. Importantly, all 
our results are confirmed numerically. 

We consider beam propagation under type-I SHG conditions in a lossless QPM 
slab waveguide. The slowly varying envelope of the fundamental wave, 
Ei = E±(x,z), and its SH, E 2 = E 2 (x,z), are (p6|) 



. dE 2 ad 2 E 2 
dz 2 dx 2 



+ %^Fr + d(z)l%eV' = 0, (2) 



where in all cases of practical interest a ~ 0.5. The normalized wave-vector 
mismatch is introduced via the real parameter (3 = kiu^Ak, where Ak = 
2ki — k 2 , and ujq is the beam width, and k\ t2 are the linear wave numbers of 
the fundamental and SH, respectively. The scaled transverse coordinate, x, 
is measured in units of ujq and the propagation coordinate, z, is measured in 
units of 21 d where Id = kiu 2 /2 is the diffraction length of the fundamental 
wave. The spatial periodic modulation of the nonlinearity is described by the 
grating function d(z) whose amplitude is normalized to 1, and whose domain 
length we define as A = tt/k, where k is the spatial grating frequency or the 
QPM frequency, which we assume real and positive. In the case of second- 
harmonic generation in lithium niobate pumped at A w ~ 1.5 /im with a beam 
width of ujq ~ 20 /mi, the intrinsic material wave vector mismatch is of the 
order \/3\ ~ 10 3 . Thus the QPM frequency must also be of this order which 
corresponds to a domain length of ~ 10 /im. For such values, a scaled grating 
frequency of the order k ~ 10 correspond to a domain length of about 1 mm. 

The two-period grating function, d(z), consists of a primary grating, d^(z), 
and a superimposed secondary grating, d^^z). We expand d(z) in a Fourier 
series 

d(z) = ^2 exp(i/cKi2;) x df^ exp(ilK 2 z) (3) 
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Fig. 2. Peak splitting in the two-period QPM grating. Dashed peaks indicate the 
location of the l'st order peaks in the one-period case. 

where the summations are over all (k, I) from — oo to oo. If we assume the 
grating functions to be square, only the odd harmonics enters into the ex- 
pansion, d 2 i+i = 2/m{2l + 1) and d 2 \ = 0. In the case of low-depletion SHG, 
the effect of the superimposed period is to split each peak of the original 
one-period QPM grating into an infinite family of peaks fll^) . More formally, 
one has peaks at all spatial QPM frequencies mux + nn 2 where m and n are 
the QPM orders related to the primary and secondary grating, respectively. 
These peaks are not delta-like, though for simplicity they are depicted as such 
in Fig. 2, but rather sinc-like functions around the relevant QPM frequency. 
The SHG-efficiency is higher the lower the order. Hence, in one-period QPM 
one would choose to match to the first peak (m = ±1). In two-period QPM 
there is an equivalent rule, i.e. the low-depletion SHG-efficiency is highest for 
m = n = ±1. To be able to treat the problem mathematically we have to 
neglect all overlap between peaks in order to avoid resonances. Thus we must 
assume k » 1 in the one-period case. In the general two-period case the 
spectrum becomes dense, as in the case of a Fibonacci grating(^). In order 
to be able to neglect any overlap between the peak we are looking at, and the 
rest of the dense spectrum, we have to assume, not only high spatial QPM 
frequencies, but also that the two frequencies are of different order. We remark 
that soliton excitation is expected to be possible no matter what the spectrum 
looks like, as long as the power is high enough fl2"T|). Our numerical simulations 
presented here agree with this expectation. It is the lack of analytical tools 
which necessitates the assumption of well separated peaks. 

By applying the asymptotic expansion (p7[) technique we have established a 
perturbation theory describing the propagation of the averaged fields in this 
two-period QPM system. We make the transformation E%{x, z) = w(x, z) and 
E 2 (x, z) = v(x, z) exp (iez), where e = j3 — mK\ — nn 2 , is the residual phase 
mismatch, which is assumed to be small. The functions w(x, z) and v(x, z) are 
assumed to vary slowly on the scale given by the QPM periods and can be 
expanded in the Fourier series 

w(x, z) = w (x, z) + J2 Mx, z V kKlZ + E ^> z ) ellK2Z i ( 4 ) 
v{x, z) = v Q {x, z) + J2 Vk{x, z)e lk ^ z + £ u t (x, z)e a ^ z . (5) 

We remark that this assumption is valid also when the QPM frequencies share 
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common harmonics. In this case the SHG spectrum discussed above becomes 
discrete as in the one-period case. However, to avoid overlap between peaks 
the QPM frequencies must still be assumed of different order. Furthermore we 
formally assume that f3 » 1, Kx » 1, and k 2 >> 1. Substituting (4-5) into 
(1-2) and matching leading order terms yields the relations for the harmonic 
coemcients([TJ). One gets 

W 1 = ^T d n )d{ (q+m) W >^ V n = — cfit^ m) ™0> (6) 

qui qt^x 

Uq = W * U0) y n = JLd^d^ W 2 Q , (7) 

qK>2 q&2 
and the equations governing the averaged fields are of the form 



1- o^TT" + VWqVo + 7(Fo| - | Vol )w = 0, 



.dw 1 d 2 w 

dz 2 dx 2 
.dv ld 2 v 

5 1 

dz 4 dx 2 



ev + rjwl - 2^\w \ 2 v = 0. (9) 



System (8-9) conserves the same power as system (1-2), P = /(|w;o| 2 -|-|fo| 2 )<ix = 
/(|-Ei| 2 + \E 2 \ 2 )dx. The strengths of the averaged nonlinearities are given by 

V = dm<$? ~ Yl d[ ^ dll_ m)! n +n \, (10) 
k 

<2) d {2) d (1) d {1) m d {1) d {2) d {2) 

In these expressions, q, I, and k are integers, and the summations are every- 
where over all integers from — oo to oo except zero. We emphasize that the 
averaged model (8-9) is valid for arbitrary values of Kx and k 2 as long as the 
assumptions Kx >> 1, and k 2 » 1 are not violated. However care must be 
taken when calculating the sums in (10) and (11). When necessary, k and I, 
themselves integers, must always be chosen such that the (k — ra)Kxl k 2 and 
(/ — n)/t 2 //ti are integers. This ensures that the Fourier coefficients d^ and d[^ 
exist. We remark that if one of the QPM frequencies is an irrational number 
expressions (10) and (11) are greatly simplified 

i7 = dSM°, (12) 

_ d^d^j y d\^ +q] ^ d$d£j y dfnlg] , 13 j 
Kx g q k 2 q q 

To make the cubic self- and cross-phase modulation symmetric in (8-9) we 
have had to assume that dffi = d_ n = —d$, which is the case for square 
grating functions. Otherwise expressions (10-13) hold for any grating function 
of the form (3). Rewriting (10) and (11) for square grating functions we get 
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E 7 x h ( 14 ) 

s {2s + m) (2s^ + n) ) 

1 = J1 (j-y I 

7T 4 2g(2r + n) (2q - 2r^ + m) (m + 2q) 

1 v.. ,. ' V («) 



™k 2 ^ 2g(2s + m) (2g - 2s*± + nj (n + 2q) 

where the summations now are over all integers q, r, and s from — oo to oo 
except zero. Notice the different sign of the lowest-order contribution to (14) 
at the (m — l,n — l)-peak, and the (m = l,n = — l)-peak. Such change 
of sign has implications to the soliton features, as shall be discussed below. 
Equations (8-9) were first derived in (|i~9|). They have the same form regardless 
of the specific type of grating, the parameters still merely being given as sums 
over the Fourier coefficients of the grating. Thus they were derived in ([24]) for 
QPM with both a linear and a nonlinear grating with or without a dc-value 
of the nonlinear grating. They were also derived in (|23| ) for a more exotic 
two-period grating. 

By letting w (x, z) = w(x,z)/r) and Vq(x, z) — v(x, z)/t) we can reduce system 
(8-9) to 

,dw ld 2 w # l2 . i2n _ 

* Tz + 2 dx^ + + l{M ~ ^ ^ = ' ( } 

Ov 1 d^v 

+ - ev + w 2 -2^\w\ 2 v = (17) 

oz 4 ox 2 

The stationary bright solitary families of system (16-17) are of the form 
w(x,z) = Ui(x) exp(zAz), v(x,z) = u 2 (x) exp(i2Xz). They are parametrized 
by the soliton parameter A > max{0, — e/2}. The stationary solutions of (8- 
9) are parametrized by exactly the same parameter but since this system is 
characterized by three (e, 7], and 7) rather than the two parameters of (16-17) 
(e and 7 = ^/rj 1 ) each solution of (16-17) corresponds to an infinite family 
of solutions of (8-9), i.e. an infinite number of QPM frequency combinations 
and thus physical setups. System (16-17) conserves the renormalized power 
P = J(\w\ 2 + |w| 2 )<ix = rj 2 P. We note that also the assumed slow variation 
of the functions w(x,z) and v(x,z) imposes constraints on the choice of the 
soliton parameter A. The characteristic soliton scale length, 2tt/X, must be 
much larger than the QPM period A = tt/k, where k in the two-period case 
is the smallest QPM frequency. 

Merely looking at expressions (14-15) suggests that in principle there is ample 
room for engineering of the averaged nonlinearities. The main parameters are 
the QPM frequencies k% and In general we can state that the higher the 
QPM frequencies are, the lower the induced averaged cubic terms become. 
Thus 7 — > for Ki,K2 — » 00. One also notices that the ratio between K\ 
and k 2 and vice versa plays an important role in the sums in expressions 
(14) and (15). We observe that the smaller the ratio the more significant 
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Fig. 3. Fraction of power in the SH for the two-period QPM system with 
(«i, k 2 ) = (997, 13) as a function of total power. Solid line and •: Theoretical and nu- 
merically measured values, respectively. Dashed curves: Zero-order approximation, 
7 = 0. The value of e is indicated at each pair of curves. 

the corresponding sums. From (14) it is obvious that the induced part of 
the averaged quadratic nonlinearity depends only on the ratio k^/ki. Thus we 
can realize the same averaged quadratic nonlinearity at large QPM frequencies 
(no induced averaged cubic terms) and at low QPM frequencies (significant 
induced averaged cubic terms). We note that if we let the second grating 
function dS 2 '(z) = 1, the only Fourier coefficient left in the expansion of S 2 '(z) 

is d^ = dffi = dll = 1 and hence 77 = d$ and 7 = ^ d^ +q Jq- The whole 
system thus degenerates into the one-period case and the results of (|TDD are 
reproduced. 

Depending on our objectives we can make the individual terms of the averaged 
nonlinearities more or less significant. Thus we can essentially make the two- 
period QPM system behave like the one-period system with induced averaged 
cubic nonlinearities by choosing K\ » k 2 while still assuming k 2 >> 1. In 
this way one can use K\ to reduce a large intrinsic phase mismatch, f3, to 
a small effective mismatch which can then be manipulated further with k 2 . 
Thus, if we let K\ — > 00 the averaged nonlinearities of for the system with 
square grating functions (14-15) with m = n = 1 reduces to 

V = ~- 2 , 7 = A( 1 -^) ( 18 ) 

7T Z K 2 7T Z \ K J 

To verify the model in this K\ — > 00 limit we have made a series of simulations 
launching soliton initial conditions, found from (16-17), in a real two-period 
structure described by (1-2) with K\ chosen to be a high prime and n 2 a 
small prime. Choosing primes effectively eliminates contributions from the 
remaining sums in (14-15) because r and s must be chosen such that tk 2 /ki 
and SK\j k 2 are positive integers. We measure the fraction of power in the SH 
after any initial transient has died out and compare this to the theoretically 
predicted value. In Fig. 3 we have summarized the results of this investigation 
for different powers and different effective phase mismatches. 

If both QPM frequencies are of the same order, the full non-simplified expres- 
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sions (14-15) must be used to calculate the averaged nonlinearities. Because 
of the way the ratio K\/k 2 enters in (14-15), one has discrete resonances be- 
tween the QPM frequencies. At these discrete resonances the induced averaged 
nonlinearities become stronger. The strength depends on the order of the reso- 
nance, the lower the order the higher the strength. To show how the strength of 
the averaged nonlinearities are calculated we concentrate on the summations 
involving sk\/k 2 . The contribution of the terms involving tk 2 /ki is handled in 
an analogous way. First-order resonances arise when k±/k 2 is a positive integer; 
then s is summed over all integers except zero from — oo to oo. Second-order 
resonances arise when K\ = (2p — 1)k 2 /2 where p is a positive integer. In 
this case s only spans the even integers. Third-order resonances occur when 
K\ = p * k 2 /3 where p — 1,2, 4, 5, 7, 8, • • •; then s must be summed over the 
integers ±3, ±6, ±9, • • •. And so on and so forth for higher-order resonances. 

We emphasize that resonances in principle must be taken into account when- 
ever the gratings periods are rational values. In actual practice, the values 
of the grating periods are not integers but may be approximated by rational 
numbers. For example, let (k,i,k 2 ) = (194.7,13.3). In this particular case, we 
can write K\ = pn 2 / 133, where p is any integer which has no common divi- 
sor with 133 except 1, showing that we are thus dealing with a 133'rd order 
resonance. In this particular case, p = 1947, and s must be summed over 
±133,±266,±399,---. 

In Fig. 4 we show the induced averaged cubic nonlinearity of the reduced 
system, 7 = ^/rf, as a function of the QPM frequency k x (k 2 is fixed). We 
have plotted the values for the lowest order resonances together with the no- 
resonance curve, i.e. the limiting curve where the average nonlinearities are 
given by (12-13) which in the case of square grating functions reduces to 



It is evident from Fig. 4 that the higher order resonances quickly converge 
towards the no-resonance curve with K\ and that only the first order reso- 
nances yield induced averaged cubic nonlinearity significantly higher than the 
no-resonance case. We note that the resonance peaks shown are truly dis- 
crete. Furthermore, when the grating periods are non-integer numbers, only 
very high-order resonances occur. The point labelled o in Fig. 4 illustrates the 
idea. It corresponds to (ki,k 2 ) = (8.66,13). While the very similar structure 
corresponding to k 2 ) = (26/3, 13) yields a third-order resonance contribu- 
tion, («!, k 2 ) = (8.66, 13) yields a 650th order resonance. Hence, the negligible 
deviation relative to the no-resonance curve. The conclusion is thus that only 
the no- resonance curve is of practical experimental interest. 

However, to verify that the model is correct we still need to analyze the full 
resonant structure of the averaged nonlinearities. Thus in Fig. 5 we focus 
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Fig. 4. Averaged cubic nonlinearity at the (m = n = l)-peak in the reduced system 
as a function of the QPM frequency k\ (k2 = 13). Full curve: No-resonance limit. ♦: 
l'st order resonances. ■: 2'nd order resonances. •: 3'rd order resonances. The lines 
between •, ♦ and a's are only to help the eye. Dashed line: K2 — > 00 approximation, o: 
Indicates a point, («i, 7) = (8.66, 0.090) very close to the p = 2 3'rd order resonance 
point («i,7) = (26/3,0.157). *: Indicates the point where 7 ps 0.050 (see text). 




Fig. 5. Averaged nonlinearities at the first order resonances from Fig. 4. K2 = 13 and 
K\ = p * 13 (p is an integer). ■: The effective cubic nonlinearity, 7, of the averaged 
system. •: Amplitude of the averaged quadratic nonlinearity, 77. The curves are 
discrete and the lines in between the «'s and a's are just to help the eye. The 
dashed and dotted lines indicate the asymptotic values of 7 (rei — > 00) and |ry|, 
respectively. indicates the point where 7 « 0.050 (see text). 

on the first order resonances and show how the effective averaged quadratic 
nonlinearity changes with the ratio k\/k2- In the same figure we have also 
plotted 7 to emphasize the interplay between 7 and rj. Recall that the graphs 
are not continuous but discrete functions of K4//C2, i-e. the plot only holds 
for values of K\ which are integers of k 2 («2 = 13 in this case). In Fig. 6 we 
have made a series of numerical experiments to confirm that the averaged 
model derived in this paper does predict correct behavior also at the QPM 
frequency resonances. Again we launch soliton found from (16-17) in the real 
system described by (1-2), but now we keep the input power fixed, while we 
change K\ in integer steps of k 2 - Simulations where carried out in the actual 
two-period QPM structure. As before we let the evolution continue until any 
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Fig. 6. Fraction of peak-intensity in SH at the first order resonances from Fig. 5. 
K2 = 13 and k% = p * 13 (p is an integer). The input power is kept constant at 
P = 15 and e = 0. The lines represent theoretically predicted values whereas ♦, • 
and ■ are numerically measured values. Dashed line and ♦: Maximum value. Solid 
line and •: Average value. Dashed dotted line and ■: Minimum value. Notice that 
the lines actually represent the discrete integer values. They are drawn as lines just 
to help the eye. 

initial transient has died out. Then we sample the minimum, maximum, and 
averaged peak intensities and compare those with the theoretically predicted 
values. 

It is worth recalling that for every soliton solution of the reduced system (16- 
17) we have a whole family of solutions of the form (4-5). In Fig. 7 we have 
shown how one specific solution for 7 « 0.050 leads to different initial condi- 
tions for two different physical setups. The two setups here chosen corresponds 
to the point -j* indicated on both Fig. 4 and Fig. 5. We also show the evolution 
of these initial conditions in Fig. 7, i.e. the evolution of the peak intensities. 

For the induced averaged nonlinearities addressed in this paper to be of po- 
tential practical importance, they have to impact the observable soliton prop- 
erties, including their excitation conditions. To show that this seems to be 
the case, in Fig. 8 we show the behavior found for the soliton content (|38"D , 
SC, as a function of the residual phase mismatch in a two-period structure 
with QPM frequencies (/ti,/t2) = (195,13). We also include the outcome of 
the numerical experiments performed for a structure with the non-integer val- 
ues («i, k 2 ) = (194.7, 13.3), to show that the effect predicted is not restricted 
to a particularly suitable choice of these parameters. We launch a FF signal, 
with a sech-shape, and no SH seeding and calculate how much of the initial 
power, P in , is bound in the soliton which eventually forms. We propagate until 
a steady state has emerged and then we collect the energy flow in a window 
wide enough to enclose substantially all the soliton. Typical values are prop- 
agation until z = 10 3 and collection of energy flow in a window of x = ±10, 
but these values were adapted whenever needed in order to capture always 
all the soliton energy. Simulations where carried out in the actual two-period 
QPM structure. The bandwidth of the SC for sech inputs can be estimated 
by using the Zakharov-Shabat scattering equations associated with the (1+1)- 
dimensional nonlinear Schrodinger equation, NLSE. With cubic terms, i.e. for 



10 



6 I 6 




5 10 5 10 

Propagation distance, z Propagation distance, z 




Transverse coordinate, x Transverse coordinate, x 



Fig. 7. Evolution of solitons in the real two-period system (1-2). Top figures 
show peak intensities. Bottom figures compare intensity profiles for the FF's (left) 
and SH's (right). Top left figure and dashed curves in bottom figures are for 
(ki,K2) = (195,13). Top right figure and solid curves in bottom figures are for 
(k±,K2) = (32.9,13)). The two initial conditions stem from the same soliton of 
system (16-17) with 7 « 0.050, e = 0, and A = 0.4. 
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Fig. 8. Soliton content for sech input FF as a function of the residual phase mis- 
match e = (3 — Ki — K2 {m = n = 1). Input power is P{ n = 50. Left plot is for 
( Ki ,k 2 ) = (195, 13) (77 ~ -0.38, 7 ~ 0.050) and right plot for k 2 ) = (194.7, 13.3) 
(77 ~ -4/vt 2 ,7 ~ 0.038). •': SC for two-period QPM. SC for purely quadratic 
model. The discrete points are the outcome of numerical experiments; the full lines 
in between are only to help the eye. The dashed and dotted lines are estimates from 
the limiting NLSE for a pure quadratic model and for a model with cubic terms, 
respectively. 

system (8-9), one gets the estimate 



ScJ^i^-fSf). (21) 

V Pin V V T + l e J 

In the plot we also include the actual value of SC that is numerically found in 
the corresponding homogeneous pure quadratic case. Compared to the pure 
quadratic model, the bandwidth for the two-period structure is found to get 
wider. This effect can be understood by comparing the approximate theoretical 
expressions for the two cases given by (21). There is also an expected significant 
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increase in the efficiency owing to the fact that in systems with induced cubic 
terms the exited solitons are closer to the initial sech-shape than in purely 
quadratic systems. To explain the "dip" in the SC in the case of (ki, k 2 ) = 
(194.7, 13.3) again one must turn to the SHG spectrum. We find that the 
(m = = 31)-peak is located exactly at e = 9.6, which is where the 

bottom of the "dip" seems to be. Thus we conclude that the efficiency of soliton 
excitation is diminished because of resonances between the (m = n = 1)- and 
the (m = —l,7i = 31)-peak. 

In Fig. 9 we plot the soliton content for the (fti,/^) = (195,13) case, but 
now we scan mismatches not only around the (m — n — l)-peak but also 
around the (m — l,n = — l)-peak. In the absence of resonances and average 
cubic nonlinearities, i.e. rj = ±4/n 2 , no differences are found in the behavior of 
soliton excitation at each peak. However, such is not the case when the soliton 
generation efficiencies can interfere with each other, as shown in Fig. 9. For 
example, in the particular case shown, one observes that soliton generation 
around the (m — l,n — — l)-peak takes place within a narrower band of 
mismatches and is less efficient than around the (m — l,n — l)-peak. This is 
because the average nonlinearities (14-15) are nonlinear functions of m and n 
and hence they change their relative strengths at the two peaks. A different 
behavior than that displayed in Fig. 9 might be obtained with different values 
of the two-grating periods and input light conditions, as shown in Fig. 10. 
When resonances are involved we note that at the (m — l,n — — l)-peak the 
induced part of the average second order nonlinearity is stronger than the no- 
resonance value 4/-7T 2 whereas it is weaker at the (m = l,n= l)-peak. Notice 
also that the soliton content vanishes in the intermediate region between the 
(m = l,n = ±l)-peaks. This is because in the intermediate region between 
both peaks, e.g., at the mismatch that would be exactly compensated with 
a single-period QPM structure, the interference between the two peaks acts 
as a periodic amplitude modulation of the effective nonlinear coefficient, with 
a period comparable with the characteristic soliton length. Such modulation 
constitutes a very strong perturbation that prevents soliton formation. 

We remark that one finds soliton generation around other peaks than the 
(m — l,n — ±l)-peak. For example, in the («i,/C2) — (195,13) case we 
observe a band (not shown in the plot) around e = 39 corresponding to the 
(m = l,7i = 3)-peak. However, soliton formation in this band is much less 
efficient than in the (m — l,n — ±l)-band and therefore we did not include 
it the plot. Similarly, other bands with higher QPM order also exist but one 
has to launch correspondingly high powers to excite solitons. 

In conclusion, we derived an averaged model for the soliton propagation in 
two-period QPM systems of the form (8-9), and we showed that the model 
gives an accurate description of stationary soliton propagation under a va- 
riety of input light and material conditions. We have shown that one can 
use one grating period of the QPM structure to reduce the intrinsic phase 
mismatch and the other grating period to induce averaged cubic nonlineari- 
ties. We found the induced averaged cubic nonlinearities to be strongest when 
one QPM frequency is a multiple of the other, but we showed that those are 
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Fig. 9. Soliton content for sech input FF as a function of intrinsic phase mismatch 
/3 with (k\,K2) = (195,13). Input power is Pj n = 50. The discrete points are the 
outcome of numerical experiments; the full lines in between are only to help the 
eye. The dashed curves are estimates from the limiting NLSE. The vertical lines 
located at f3 = 182 and (5 = 208 indicate the m = 1, n = —1 and m = n = 1 peaks, 
respectively. The dashed vertical line at (3 = 195 indicates the location of the peak 
in the absence of a second period (see Fig. 2). 
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Fig. 10. Same as in Fig. 9, but here (ki,«2) = (195,39) 



strictly discrete resonances that cannot occur when experimentally feasible 
grating periods are taking into account. However, our investigations predict 
that also the no-resonance averaged contributions can be important. In par- 
ticular, we found that the QPM engineered averaged cubic nonlinearities, in- 
duced in feasible two-period samples, enhances the peak-efficiency and residual 
mismatch-bandwidth of the soliton excitation process with non-soliton single 
frequency pump light. 
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